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We study thermodynamic and dynamic properties of a rigid model of the fragile glass forming 
liquid orthoterphenyl. This model, introduced by Lewis and Wahnstrom in 1993, collapses each 
phenyl ring to a single interaction site; the intermolecular site-site interactions are described by 
the LJ potential whose parameters have been selected to reproduce some bulk properties of the 
orthoterphenyl molecule. A system of A'^ = 343 molecules is considered in a wide range of densities 
and temperatures, reaching simulation times up to 1 /iS. Such long trajectories allow us to equilibrate 
the system at temperatures below the mode coupling temperature Tc at which the diffusion constant 
reaches values of order 10~^° cm? / s and thereby to sample in a significant way the potential energy 
landscape in the entire temperature range. Working within the inherent structures thermodynamic 
formalism, we present results for the temperature and density dependence of the number, depth and 
shape of the basins of the potential energy surface. We evaluate the total entropy of the system by 
thermodynamic integration from the ideal -non interacting- gas state and the vibrational entropy 
approximating the basin free energy with the free energy of 6 A — 3 harmonic oscillators. We evaluate 
the configurational part of the entropy as a difference between these two contributions. We study 
the connection between thermodynamical and dynamical properties of the system. We confirm that 
the temperature dependence of the configurational entropy and of the diffusion constant, as well as 
the inverse of the characteristic structural relaxation time, are strongly connected in supercooled 
states; we demonstrate that this connection is well represented by the Adam-Gibbs relation, stating 
a linear relation between log-D and the quantity 1/TSc- This relation is found to hold both above 
and below the critical temperature Tc — as previously found in the case of silica — supporting the 
hypothesis that a connection exists between the number of basins and the connectivity properties 
of the potential energy surface. 

PACS number(s): 



I. INTRODUCTION 

Understanding the dynamic and thermodynamic prop- 
erties of supercooled liquids is one of the more challeng- 
ing tasks of condensed matter physics (for recent re- 
views see Refs. and references therein). A signif- 
icant amount of experimental |^|-^, numerical pO[ | and 
theoretical work |pl]-p^ is being currently devoted to 
the understanding of the physics of the glass transition 
and to the associated slowing down of the dynamics. 
Among the theoretical approaches, an important role has 
been played by the mode coupHng theory (MCT) |]ll| , |l2| , 
which, interpreting the glass transition as a purely dy- 
namical phenomenon, has constituted a significant tool 
for the interpretation of both experimental |p|,p|,p^|"[l9|] 
and numerical simulation results |po|-^^ in weakly su- 
percooled states. 

In the last years the study of the topological struc- 
ture of the potential energy (hyper-) surface (PES) p^ ] 
and the connection between the properties of the PES 
and the dynamical behavior of glass forming liquids has 
become an active field of research. Building on the in- 
herent structure (IS) thermodynamic formalism proposed 



long time ago by Stillinger and Weber |^ , the PES can 
be uniquely partitioned in local basins and properties of 
the basins explored in supercooled states (average basin 
depth and basin volume), have been quantified. Stud- 
ies have mainly focused on two fundamental questions: 
(i) which are the basins relevant for the thermodynamics 
of the system, i.e., which are the basins populated with 
largest probability? and (ii) which arc the topological 
properties of the regions of the PES actually explored 
by the system during its dynamics? From this point of 
view, the PES approach has somehow unified, at least on 
a phenomenological level, the thermodynamic and dy- 
namic approaches to the glass transition. 

Numerical analysis of the PES has shown that tra- 
jectories in configuration space can be separated into 
intra-basin and inter-basin components [ p^p5| . The time 
scales of the two components become increasingly sepa- 
rated on cooling. The intra-basin motion has been as- 
sociated with the high-frequency vibrational dynamics, 
while the structural relaxation (a-relaxation) has been 
related to the process of exploration of different basins. 
It has also been shown that on lowering T, the system 
populates basins of lower and lower energy p6[ . The T 
dependence of the depth of the typical sampled basins fol- 
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lows a 1/T law [p7|^9[ for fragile liquids, and, for strong 
liquids, it appears to approach a constant value on cool- 
ing |3^. The number of basins O as a function of the 
basin depth e/5 has also been recently evaluated for a 
few models [|l],|2|j28|j3^||,|^ , opening the possibility of 
calculating the so-called configurational entropy Sc and 
its T dependence. Sc, defined as the logarithm of the 
number of accessible basins Sc = ks log 51, has been suc- 
cessfully compared with theoretical predictions p3| , ^ . 
At the same time, the approaches and the techniques de- 
veloped for the analysis of the PES of structural glasses 
have spread to the field of disordered spin systems, where 
similar calculations have been performed ||35[] and simi- 
lar conclusions have been reached. The evaluation of Sc 
for models of glass-forming liquids allows to numerically 
check, in a very consistent way, the relation between Sc 
and the systems characteristic time r, proposed by Adam 
and Gibbs Q , and recently "derived" in a novel way Q . 
Numerical support for a relation between the T depen- 
dence of Sc and the T dependence of r, although limited 
to very few models, is providing new physical insight on 
the connection between thermodynamics and long time 
dynamical properties. The ideas developed within the 
inherent structure formalism have also been generalized 
to out-of-equilibrium conditions where the slow aging dy- 
namics has been interpreted as the process of searching 
for basins of increasingly deep energy ]3^^-|40|] . 

In this paper we study the properties of the PES for 
a rigid model (LW) of the fragile glass former orthoter- 
phenyl, first introduced by Lewis and Wahnstrom |4l| ] 
and recently revisited by Rinaldi et al. We have 

studied the properties of the PES in a temperature range 
in which the diffusion coefficient varies by more than 
four orders of magnitudes for five different density values. 
This work attempts to build a bridge between models of 
more direct theoretical interest, like Lennard Jones (LJ) 
and soft spheres, and models which appear to reproduce, 
even if in a crude way, properties of complex materials. 
In this respect, orthoterphenyl is the best candidate, be- 
ing one of the most studied glass forming liquids [pT| . 
The LW model is a three-sites model, with intermolec- 
ular site-site interactions described by the LJ potential. 
This model is among the simplest models for a nonlinear 
molecule. The limitation constituted by the fact that it 
does not take into account the internal molecular degrees 
of freedom (see for a more realistic model), is over- 
ruled by the observation that its simplicity — it can be 
considered as an atomic LJ with constraints — allows one 
to reach simulation times of the order of ^ls. Hence a sig- 
nificant sampling of the PES in a large temperature and 
density range is possible. Moreover, this model consti- 
tutes an ideal bridge between simple atomic models and 
molecular models, being possible to treat it under several 
approximations p^ . 

The paper is structured as follows: In Sec. ||we briefly 
recall the main results of the IS formalism. In Sec. Ill 



we show the calculation of the configurational entropy 
as a difference between the total entropy and the vibra- 
tional entropy. In Sec. IV we give some numerical details. 
We present our results in Sec. 0, which is divided into 
subsections detailing the calculation of the total entropy 
by thermodynamical integration from the ideal gas state, 
the study of the vibrational properties of the PES and 
the calculation of the configurational entropy. In the end 
we study the link between configurational entropy and 
the diffusion constant, investigating the validity of the 
Adam-Gibbs equation. In Sec. ^ we finally discuss our 
results and we draw some conclusions. In Appendix A 
we report the analytical calculation of the total entropy 
of a system of LW molecules in the non-interacting "ideal 
gas" limit. 



II. INHERENT STRUCTURE 
THERMODYNAMICS FORMALISM 

In this section we briefly review the IS formalism in 
the NVT ensemble [||,|4| , the extension to the NPT en- 
semble poses no particular problems |23| . This formalism 
has become an important tool in the numerical analysis of 
classical models since it is numerically possible to calcu- 
late in a very precise way the inherent structures (deflned 
as the local minima of the PES) and hence compare the 
theoretical predictions with the numerical results. Given 
an instantaneous conflguration of the system, a steep- 
est descent path along the potential energy hypersurface 
defines the closest IS. 

In the IS formalism, the partition function of a sys- 
tem is written as a sum over all the PES basins. Basins 
of given IS energy contribute non-negligibly to the total 
sum if their IS energy is very low, if their volume is very 
large and/or if they are highly degenerate, i.e. several 
basins are characterized by this IS energy. This corre- 
sponds to partition the phase space in the local energy 
minima of the PES and their basins of attraction. Such 
a partition is motivated by the fact that in supercooled 
states, the typical time scales of the intra basin and inter 
basin dynamics differ by several orders of magnitude and 
hence the separation of intrabasin and interbasin vari- 
ables becomes meaningful. 

In the 6iV-dimensional configuration space, the parti- 
tion function Z for a system of N rigid molecules can be 
written as: 



Z 



X3N 



J dq~exp(-I/(q^)/fcBT), 



(1) 



where denotes the positions and orientations of 
the molecules, V^(q^) is the potential energy , 7^, 
with /i = X, y, z, are the principal moments of iner- 
tia of the molecule, A^ = (27r/^fcBT)^/^//i, and A = 
/i(27rTOfcsT)~^/^ is the de Broglie wavelength. 
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Let rt{Ejs) denote the number of minima with energy 
Ejs, and f{T,Eis) the average free energy of a basin 
with basin depth Ejs- f{T,Ejs), which takes into ac- 
count both the kinetic energy of the system and the local 
structure of the basin with energy Ejs , is defined by: 



f{T,E 



IS) 



1 



A3^^ n{Eis) 



E 

basins 



dq^ exp{~{V-Eis)/kBT) 



where the configurational entropy Sc{Eis) has been de- 
fined as 



S,{Ejs) = kB ln[n{Eis)]. 
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df 
dE 



T 



IS 



dE 



0. 



IS 



eis{T) cx 



T 



III. EVALUATION OF THE 
CONFIGURATIONAL ENTROPY 

The Eq. provides a natural starting point for a 
numerical evaluation of the configurational entropy. In- 
deed, the free energy F{T, V) per molecule can be split 
in the usual way as a sum of an energy and an cntropic 
contribution. Considering Eq. (|^) we write: 



(2) 



where Rbasin is the configuration volume associated with 
the specific basin. The partition function can then be 
rewritten as a sum over all basins in configurational 
space, i.e. 

n{Eis) exp ( — 1 



V- ( -TSc{Eis) + Eis + fiT,Eis) \ 

His 



(4) 



In the thermodynamic limit, the free energy of the liq- 
uid can be calculated using 

F [eis{T)] = eis{T) + f [T, eis{T)] - TS^ [eisiT)] , (5) 

where eis(T), the average value of the IS energy at tem- 
perature T, is the solution of the saddle point equation 



(6) 



The liquid free energy expression Eq. (||) has a clear in- 
terpretation. The first term in Eq. (|^) takes into account 
the average energy of the PES minimum visited, the sec- 
ond term describes the volume of the corresponding basin 
of attraction and the kinetic energy, and the third term 
is a measure of the multiplicity of the basin. 

It can be rigorously shown ||2^,^,^ that, if the den- 
sity of state il{Eis) is Gaussian, and if the basins have 
approximately the same shape or are, to a good degree, 
harmonic, the important relation holds 



(7) 



On lowering T, basins with lower Ejs energies and lower 
degeneracy are populated, i.e., both e/5 and Sc decrease 
with T. 



F{T) = E{T) - TS{T) 

= -TSaiT) + eis{T) + E,{T) - TS,{T) 



(8) 



where the index v indicates the vibrational quantities 
(intra-basin components). In order to evaluate these 
quantities we calculate the basin free energy as the free 
energy of 6A^ — 3 independent harmonic oscillators [3^ ] 
plus a contribution that takes into account the basin an- 
harmonicities. Then we can write 

EiT) ^U^^\^ + eis{T) + C/anh(r), (9) 



5(T) = 5„(r) + 5,(T) 

= 5harm(T)-f 5anh(T)-f5,(T), 



and 



•5*11 arm 



h UJn(T) 



(10) 



(11) 



where the frequencies LOn are the square root of the eigen- 
values of the Hessian matrix calculated in the inherent 
structures. 

Thus, the total entropy is the sum of two contribu- 
tions: Sc{T) which accounts for the multiplicity of basins 
of depth eis{T), and Sy{T) which accounts for the "vol- 
ume" of the basins. The last equations give us, in a very 
transparent way, the physical meaning of the partition of 
the PES; moreover, they provide us a very efficient way to 
calculate the configurational entropy as a difference be- 
tween the total energy of the system and the vibrational 
entropy. 

The total entropy S can be evaluated via thermody- 
namic integration, starting from a known reference point. 
Every variation of total entropy can be generally written 
as the sum of variation along isochores and isotherms in 
the form: 



AS* = AS- 



V 



AS7 



(12) 



Then the change of entropy along an isochore between 
two temperatures T and T is 



ASv = S{V, T) - Siy, T) 



(13) 



and the change along an isotherm between two volumes 
V and V is 
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ASt = S{V,T)-S{V,T) 

1 



(14) 



T 



E{V,T)- E{V,T) 



dVP(y,T) 



In the present case, to evaluate the total entropy of the 
liquid we start from the known expression of the ideal 
gas of LW molecules, reviewed in Appendix A. To evalu- 
ate the basin free energy / [T, eis{T)], we select as refer- 
ence point the free energy of {6N — 3) independent har- 
monic oscillators (whose distribution of frequencies can 
be calculated evaluating the eigenvalues of the Hessian 
matrix evaluated in the IS structure) and add corrections 
to take into account the basin anharmonicities.The har- 
monic contribution to the entropy is given by Eq. (^ij). 

Assuming that the anharmonic contribution is inde- 
pendent from the basin depth, the anharmonic correc- 
tions to the entropy at T can be calculated integrating 
the quantity dUi^nh/T, where C/anh is implicitly defined in 
Eq. (|), from T = to T (sec Eq. (Ilj 



second set of simulations at constant T [T = 5000 K) in 
the volume range 10^ — 10^ nm'^ is performed to calculate 
the excess pressure ( i.e. the pressure beyond the ideal 
gas contribution). 

To calculate the inherent structures visited in equi- 
librium we perform conjugate gradient minimizations to 
locate the closest local minima on the PES. We use a 
tolerance of 10~^^ kJ/mol in the total energy for the 
minimization. For each thermodynamical point we min- 
imize at least 100 configurations and we diagonalize the 
Hessian matrix of at least 50 configurations to calculate 
the density of states. The Hessian is calculated choosing 
for each molecule the center of mass and the angles asso- 
ciated with rotations around the three principal inertia 
axis as coordinates. 



V. RESULTS 



A. Dependence of the total entropy on T and p 



IV. NUMERICAL DETAILS 



The LW model is defined by the spherical potential 



V{r) = 4e 



(7)" 



+ Ai + A2 r, 



(15) 



with e = 5.276 kJ/mol, a = 0.483 nm Ai = 0.461 kJ/mol 
and A2 = —0.313 kJ/{mol nm). The parameters of the 
potential are selected to reproduce some bulk properties 
of the OTP molecule such as the temperature de- 
pendence of the diffusion coefficient and the structure. 
The values of Ai and A2 are selected in such a way the 
potential and its first derivative are zero at Tc = 1.2616 
nm. Such a potential is characterized by a minimum at 
r = 0.542 nm of depth —4.985 kJ/mol. The integration 
time step is 0.01 ps. The shake algorithm is implemented 
to account for the molecular constraints. 

We study a (N,V,E) system composed by = 343 
molecules (1029 LJ interaction sites) at 5 different densi- 
ties (see Table |) for several temperatures at each density 
(Table ||). The total simulation time is quite long, ex- 
ceeding 10 ^s. We take care to check the thermalization 
of the system at the lowest temperatures; the production 
run follows a thermalization run whose length is compa- 
rable to the relaxation time at the considered thermo- 
dynamical point. We are able to thermalize the system 
at temperatures at which the diffusion constant reaches 
values very low, of order 10~^° cm^/s, 10^ smaller then 
at ambient temperature. 

Two additional simulations are performed to connect 
the range of densities and temperature studied with the 
ideal gas reference point. The system at density is 
simulated for temperatures ranging from 250 to 5000 K 
to evaluate the T dependence of the potential energy. A 



To estimate the total entropy for the model we proceed 
in three steps as shown in Fig. ^ The thermodynamic 
path has been chosen to avoid the liquid-gas first order 
line. 

(1) Integration along the isotherm Tq — 5000 K from 
{To,V = 00) (perfect gas) to (Tq, V4=118.35 nm^), cor- 
responding to point Co in Fig. |^. The ideal gas contri- 
bution to the total entropy is discussed in Appendix A. 
The entropy at Co can be calculated as 



U(To,V4) f^^ dV 



To 



To 



(16) 



where P^^^ is the pressure that exceeds the pressure of the 
ideal gas, i.e. the contribution to the pressure due to the 
interaction potential and U is the system potential en- 
ergy. The values of the pressure Pex{T = To,V, N = 343) 
as a function of V are reported in Fig.^ (a). Pex{T — 
To,V,N — 343) has been fit using the virial expansion 



P,,(T ^To,V,N = 343) = ^ a^F-f^+i) . 



(17) 



fc=i 



The Ofc values are reported in Table III, from which we 



estimate the first virial coefficient at To 

B2{To) = ai/iksToN^) = 0.596 nm^ (18) 

In Fig.^ (b) we plot the potential energy as a function of 
volume along the T = To isotherm. 

The total entropy at the reference point Co is S{Co) — 
294.8 J/( mol K ), resulting from the sum of three con- 
tributions 
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^.^(Co) = 339.03 J/{molK), 
f^" dV 

/ —PeAV,To)^-U.9J/{molK), 

Joo ^0 



(19) densities is reported in Table [V . These values are used as 



and 



UjCo) 
To 



= 0.64 J/{molK) 



(20) 



(21) 



(2) Integration along the isochorc V ~Va from To to 
T* = 380 K, corresponding to the point C'l in Fig.|l|. To 
evaluate the entropy along this isochore we use 



SiT*,V4) = S{To, Vi) + 3R log(T7To) 



T* 



dTdU{V4,T) 



T 



dT 



(22) 



Fig. H (a) shows the potential energy for the V — Vi iso- 
chore. To calculate the integral in Eq. (|2^), we fit the 
potential energy using the functional form which best in- 
terpolates the calculated points 



U{Vi, T)=uo + ?iiT3/5 + U2T, 



(23) 



obtaining the values uq = — 94.405,mi = 0.533, U2 = 
0.00183 (energy in kJ/mol). 

The total entropy at the reference point Ci is S{Ci) — 
191.8 J/( mol K ), resulting from the sum of three con- 
tributions: 



S'(Co) = 308.6 J/{molK), 

3i? log(r/5000) = -64.3 J/{molK), 



and 



T* 



dTdU{V4,T) 



To T dT 



-52.5 J/{molK). 



(24) 
(25) 



(26) 



(3) Integration along the isotherm T* from V4 to a 
"generic" V. To determine the total entropy difference 
for all studied densities we calculate 

S{T*,V) - SiT*,V4) - Sm{T*, V) - S^T*, V^) 

1 



+ — [UiT*,V)-U{T*,Vi)] 

dV 

+ -^PeATW). (27) 



Figs. 1^ (b) and (c) show respectively the potential en- 
ergy and the excess pressure as a function of volume at 
T = T* . For convenience we fit P^x with a third order 
polynomial 



(28) 



k=l 



where the values of the coefficients are given in Ta- 
ble in. The resulting total entropy at T* for all studied 



reference entropies for the T dependence of S. For each 
of the studied isochores, we calculate the T-dependence 
of the total entropy according to Eq. (|2|). In this low 
T-range, the potential energy is very well represented by 
the Rosenfeld-Tarazona law fE5| 



U{V,T) ^ Uo{V) + a{V) T^/^ 



(29) 



consistent with what was found for LJ systems. In Fig. ^ 
we show the temperature dependence of the potential en- 
ergy at all densities. The best-fit Uq{V) and a{V) values 
are reported in Table 0. 

The calculated total entropies at each considered den- 
sity are plotted in Fig. ||. 



B. Dependence of the inherent structure energies on 
T and p 



In Fig. g we show the temperature dependence of the 
energy of the calculated inherent structures together with 
a fit (according to Eq. (0)) in the form 



eis{V,T) ^ Aiy) + 



BjV) 
T 



(30) 



The values of the fitting coefRcients A{V) and B{V) are 
reported in Table On lowering temperature the sys- 
tem populates minima of lower and lower energy. It is 
worth noting that, in contrast to the case of the actual 
potential energy, the slope of these curves varies strongly 
with densities. 

From the T and V dependence of ejs the anharmonic 
potential energy can be calculated. Fig. |^ shows t/anh(T) 
for two densities (symbols). We also show a cubic ex- 
trapolation (solid lines) in the form of 



(31) 



As shown in Fig. |^, the anharmonic contribution is rather 
small, in agreement with previous findings for the LJ 
model. For this reason, the low signal to noise level does 
not allow a well-defined characterization of the C2 and C3 
values. To decrease the number of free parameters, we 
consider C2 to be volume independent, and we fit simulta- 
neously, according to Eq. ( |3T| ) , C2 and the V dependence 
of C3. As we will show in the following, the anharmonic 
contribution to the entropy is much smaller than the har- 
monic one and hence the choice of C2 and C3 does not 
affect significantly the resulting configurational entropy 
estimate. 
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C. Density of states and vibrational harmonic 
entropy 



isotherms increases decreasing density. The hnes are in- 
terpolations of the data using the fits of Fig. |. 



In this section we study the shape of the basins by 
investigating the properties of the density of states and 
we calculate the vibrational harmonic entropy. In Figs. ^ 
(a) and (b) we show the temperature and density de- 
pendence of the density of state, namely the histogram 
of the square root of the eigenvalues of the Hessian cal- 
culated for the inherent structures. The distribution is 
characterized by only one peak, not showing any clear 
separation between translational and rotational dynam- 
ics; the width of the distribution increases on increasing 
temperature. The position of the maximum is found to 
be to a good extent independent of temperature; at vari- 
ance it increases with density as the width does. These 
features show that the LW PES basins have shapes that 
are function of the energy depth and of the density. 

It is worth noticing one particular feature of Fig. |^ (b); 
all the curves cross at a value of the frequency oj* ~ 44 
cm~^. The presence of this isosbestic frequency (in anal- 
ogy with the well-know isosbestic frequency observed in 
the Raman spectrum of water |4^ ) supports the possibil- 
ity that a two-state model |Q may provide a reasonable 
description of the change of the density of states with 
temperature and, correspondingly, of the change of the 
density of states with the basin depth. 

In Fig. ^ (a) and (b) we plot the quantity 
-^"^ -^'^s(^fc/^o) as a function of T and of the 

e/5 respectively. The scale frequency uJq is chosen as 
1 cm^^. This quantity is an indicator of the average 
curvature of the basins and, being a sum of logarithms, 
is very sensitive to the spectrum tails. As shown in 
Fig. H (a) Yl^k=i^ log(ti^fe/ti-'o) increases with temper- 
ature along isochores and increases with density along 
isotherms. 

As noted previously for the LJ |^,^ and for the 
simple-point charge extended (SPC/E) model for wa- 
ter H^, the dependence of Yl\=i^ ^og{oJk/iOo) from 
e/s can be well approximated by a linear dependence, i.e. 



1 



fe=i 



h UJn{T) 



^a{V)+b{V)eis{T), (32) 



where To defines the T scale (To = IK) . This dependence 
indicates that deeper and deeper basins have larger and 
larger volumes (their average frequency being smaller). 
The fact that basins of different depths have different 
volumes introduces an important contribution to Eq. (|^) 
since the term df /deis is different from zero. The impli- 
cation of this non-zero contribution has been discussed 
recently in Refs. |3,||,|§. 

In Fig. ^ we show the harmonic contribution to the 
entropy as calculated from Eq. (^ij). This contribu- 
tion is obviously increasing with temperature and along 



D. Vibrational anharmonic entropy 

Integration of the anharmonic energy f/anh, obtained 
from Eq. (||) according to Eq. (|l3|) , gives directly the an- 
harmonic contribution to the entropy. For the LW case, 
Uanh is described by the polynomial in T of Eq. (|l]) , and 
we obtain 



^anh(T) = 2c2T-f -CaT^. 



(33) 



The inset of Fig. shows the anharmonic contribution to 
the vibrational entropy as calculated by integrating the 
anharmonic contribution to the potential energy. This 
contribution is negative showing that, in the range of 
densities and temperatures studied, the leading anhar- 
monic contribution acts in the direction to decrease the 
volume of the basin. 



E. The configurational entropy 



In Fig. 11 we plot the configurational entropy calcu- 



lated subtracting the vibrational (sum of the harmonic 
and anharmonic terms) from the total entropy for the 5 
studied isochores. As expected the degeneracy of basins 
increases on lowering density, in agreement with the ev- 
idence that a glass transition may be induced along an 
isothermal path by progressively increasing the pressure. 
Considering Eqs. (|ol) (0) (|^) (||)(|3l), the configura- 
tional entropy can be described in the entire density and 
temperature range considered by means of the functional 
form 



Se{T) = 5(r) - 6 



a{V) + b{V) 



N 



AiV) 



B{V) 



T 



(34) 



2C2T - -cgT-^ 



These curves are plotted in Fig. as solid lines. In 
the range of temperatures and density studied, Sc/R per 
molecule varies from about 4 to 3, a figure not very dif- 
ferent from the estimated configurational entropy of or- 
thoterphenyl, based on an analysis of the T dependence 
of the measured specific heat [^^. We recall that the 
LW model represents each phenyl group as one single in- 
teraction site and it does not account for the the molecule 
fiexibility. The similar estimate of Sc seem to suggest 
that steric effects are dominant in controlling the config- 
urational entropy. 
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F. Diffusion and the Adam-Gibbs relation 

In order to investigate the connection between the long 
time dynamics of the system and the underlying PES, 
we calculate the center-of-mass diffusion coefficient D{T) 
from the mean-square displacement via the Einstein re- 
lation 



D{T) ^ lim l(r2(t)> 



(35) 



To guarantee a proper diffusive regime, at all densities 
simulations are performed until the average mean square 
displacement is greater than 0.1 nm^ at the lowest tem- 
peratures and 10 nm^ at the highest. The inverse of the 
diffusion coefficient provides an estimate of the charac- 
teristic structural relaxation time of the LW model. 

The D values calculated are shown in Fi g, Fig. |l^ 
(a) shows the dependence on T , while Fig. |13| (b) shows 
the dependence on 1/T. Fig. |l^ (a) also shows the best 
fits to the power law 



D{T) cx (T - Tc)'^ 



(36) 



predicted by the ideal MCT in weakly supercooled states. 
The consistency of the MCT prediction for a wide range 
of D values confirms the analysis of Rinaldi et al. ji^ ] 
where explicit ideal MCT calculations were presented and 
successfully compared with the numerical results along 
one isobar. Fig. |l^ shows also that clear deviations from 
the ideal MCT take place when the diffusion value be- 
comes smaller than lO^^cm^/s. The representation of 
D as a function of 1/T shown in part (b) shows that 
the ideal MCT region is followed by a T region where 
new types of processes become effective in controlling the 
molecular dynamics. These processes, termed hopping 
processes, transform the ideal MCT divergence of char- 
acteristic times into a crossover. In the region of D values 
between 10~*cm^/s and 10~^°cto^/s, limited from below 
by the present numerical resources, data are consistent 
with an apparent Arrhenius dependence with parame- 
ters which could well become T dependent if studied in 
a larger range of D values. 

The ideal MCT critical temperatures and 7 values, de- 
termined by the fit of the D values to Eq. (^6|), as a 
function of density are shown in Fig. The density 
dependence of Tc is almost linear. The exponent 7 seems 
to increase on increasing density, but the noise does not 
allow us to rule out the possibility of a constant value. 
The filled circle indicates the value of the critical tem- 
perature Tc = 265 K determined from an isobaric run in 
Ref. 

We finally study the link between configurational en- 
tropy and diffusion coefficient, investigating the validity 
of the Adam-Gibbs equation. Fig. |l5| shows log 13 as a 
function of l/{TSc); for all studied isochores, log 13 vs. 
l/{TSc) is well described by a linear relation, with coef- 
ficients which are volume dependent, as previously found 



for the LJ Hquid [|3,for the SPC/E model for water H] 
and for the BKS model for silica 

We note on passing that deviations from linear behav- 
ior are observed at large value of log D, where intra- and 
inter- basin dynamics time scales are no longer separated. 
At high T, it has been proposed [|2j that entropy — as 
opposed to configurational entropy — is the relevant ther- 
modynamic quantity controlling dynamics. 



VI. DISCUSSION AND CONCLUSIONS 

In this article we have studied systematically the prop- 
erties of the potential energy surface for a simple three- 
site rigid model designed to mimic the properties of the 
fragile glass forming liquid ortho-tcrphcnyl. The choice 
of this simple model, which collapses the entire phenyl 
ring into one interaction site, allows us to run very long 
trajectories and to study in supercooled states the molec- 
ular dynamics up to 1 /is, allowing the determination of 
diffusion coefficients down to 10~^°cm^/s. 

Wc have found that, as in the atomic LJ case, by cool- 
ing along an isochore, basins of the PES of deeper and 
deeper energy are explored. The basin volumes are func- 
tions of the depth in agreement with previous studies. 
Using the inherent structure thermodynamic formalism, 
we have calculated the number of basins of the PES and 
their depth, in the region of depth values probed by our 
simulations. As a result, we presented a full characteri- 
zation of the the temperature and density dependence of 
the basin depth, degeneracy and volumes. 

These results are used to provide a consistent model for 
the intra-basin vibrational entropy. This, together with 
the numerical calculation of the total entropy via ther- 
modynamic integration starting from the ideal gas state, 
allow us to calculate the configurational entropy — the 
difference between the total entropy and the vibrational 
one. This quantity is of primary interest both for com- 
paring with the recent theoretical calculations |l^,Q and 
both to examine some of the proposed relation between 
dynamics and thermodynamics |^ , |l4| , ^ connecting a 
purely dynamical quantity like the diffusion coefficient 
to a purely thermodynamical quantity (S'c). To examine 
such a possibility we compare for five different isochores 
the T dependence of D with the Adam-Gibbs relation. 
In the entire range of T and densities studied the Adam 
Gibbs relation appears to provide a consistent represen- 
tation of the dynamics for the LW model. 

It is important to observe that a linear relation between 
log 13 and l/{TSc) holds both above and below the ideal 
MCT critical temperature T^, in agreement with similar 
finding for the silica case ||3^ . Recent works based on the 
instantaneous normal mode technique for several rep- 
resentative models p5|~|59| provides evidence that above 
Tc the system is always located in region of the PES close 
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to the border between different basins. The number of 
diffusive directions significantly decreases above Tc and, 
if only data above Tc are considered, the number of dif- 
fusive directions would appear to vanish at Tc- Hence 
dynamics above T,. is a dynamics of "borders" between 
basins and there is no clear reason why such dynamics 
should be well described by the Adam-Gibbs relation, 
which focuses on the "number" of basins explored. The 
observed validity of the AG relation — both above and 
below Tc — reported in this manuscript supports the hy- 
pothesis that a direct relation exists between the number 
of basins and their connectivity 59 1 . It is a challenge 
for future studies to confirm or disprove this hypothesis. 
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so the total partition function for an ideal gas of OTP 
molecules can be expressed as 



Nl 



(41) 



We approximate A^! w N'^e^^ . The free energy F^^j and 
the entropy Sid of the non-interacting system then be- 
come 



F,d{T, V, N) = -fcsTln [Z,d{T, V, N)] = 



N 



- In 2 + In V"y^y^37^^7^2;7^^ + 3 In T - In Af + 1 

S,a{T, V, N) = - i^,^(T, V, N) 

Kb oi 



NkB < 4 + i In TT — ln2 + In 



A^ 



(42) 



(43) 



where the term In 2 is due to the two possible degenerate 
angular orientations of the molecule pC| . 



APPENDIX A: IDEAL GAS ENTROPY FOR 
THE LW MODEL 



In this appendix we calculate the partition function of 
a system of A'^ LW molecules in the non-interacting - 
ideal gas- case. The LW OTP molecule Q is a rigid 
three sites isosceles triangle; each site represents an en- 
tire phenyl ring of mass m = Qmc — 78 a.m.u., where 
mc is the mass of the carbon atom, and it is consid- 
ered as the LJ interaction site. The length of the two 
short sides of the triangle is cr = 0.483 nm and the angle 
between them is 9 = 57r/12 (75 degrees). 

The three moments of inertia for the single molecule 
are: 



Ix — —mcr cos 
o 



= 1.248 X 10"'''* kg m^, 



ly = 2ma^ sin^ ( ^ ) = 2-204 x 10"'"' kg (37) 



and 



t: cos , „ 
3 V2 



2 sin^ 
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k 








1 


1.036 


126.647 


5.022 


2 


1.060 


123.883 


4.985 


3 


1.083 


121.120 


4.948 


4 


1.108 


118.356 


4.910 


5 


1.135 


115.593 


4.871 



TABLE I. Densities, volumes and simulation box lengths 
calculated. 



Pi 


P2 


P3 


Pi 


P5 


170 


190 


230 


280 


320 


185 


200 


240 


300 


340 


190 


210 


260 


320 


360 


195 


230 


280 


340 


380 


210 


250 


300 


360 


400 


220 


280 


320 


380 


420 


240 


300 


340 


400 


440 


260 


320 


360 


420 


460 


280 


340 


380 


440 


480 


300 


360 


410 


460 


530 








480 





TABLE IL Temperatures (in K) for which calculations are 
performed. 



i 


tti (MPa nm^''+i)) 


Pi (MPa nm^('+i') 


1 


4835.96272 xlO'' 


15943.2 


2 


1000.53765 xlO" 


-256.591 


3 


9654.69470 xlO® 


1.1745 


4 


3873.87001 xlO^° 


-0.00111551 



TABLE III. Fitting coefficients for the excess pressure as 
a function of 1/V at T = 5000 K and at T = 380 K 



k 


S{T*) J/{molK) 


1 


192.80 


2 


188.21 


3 


183.54 


4 


177.95 


5 


172.12 



TABLE IV. Total entropy at five densities for the reference 

temperature T*. 



Pk 


Uo (kJ / mol) 


a (kJ T-^^V mol) 


A{kJ/mol) 


S (kJ T / 


1 


-86.30 


0.4385 


-79.11 


-285 


2 


-88.94 


0.4716 


-80.14 


-436 


3 


-92.07 


0.5231 


-81.88 


-676 


4 


-95.23 


0.5762 


-81.36 


-965 


5 


-96.06 


0.5731 


-81.89 


-1100 



TABLE V. First two coloumns are the coefficients for 
the potential energy U{T,V) = Uo{V) + a{V)T^^'^; second 
two columns are the coefficients for the inherent structures 
e,s{V,T) = A(V) + B{V)/T. 



pk 


a{V) 


h{V){mollkJ) 


1 


47.1 


0.342 


2 


41.2 


0.259 


3 


36.5 


0.192 


4 


32.1 


0.132 


5 


28.9 


0.869 



TABLE VI. Coefficients of the fit to the form 
YZf log('-'fcM) = <V) + b{V) eis{T). 
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FIG. 1. Thermodynamic integration paths used to calcu- 
late the total entropy at the thermodynamical points of in- 
terest starting from the ideal — non interacting- gas state. 
Details are given in the text. 
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FIG. 2. (a) Excess pressure at T = 5000 K as a function 
of volume. The open circles are the MD results. The dashed 

line is the the first term of the virial expansion to the excess 
pressure; the solid line is a third order polynomial fit to the 
entire sot of data, (b) Potential energy at T = 5000 K as a 
function of volume. 
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FIG. 3. (a) Integration Step 2. Potential energy (open 
circles) at the density p4 in the entire temperature range con- 
sidered; the solid line is the fit of the data to Eq. (23). The 
inset shows the lowest temperatures region in order to stress 
the accuracy of the fit. (b) and (c) Integration Step 3. Po- 
tential energy (b) and pressure (c). 
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FIG. 4. Potential energies at the different densities as a 
function of T^/^. The straight sohd Unes show the validity of 
the Rosenfeld-Tarazona law Eq. (29). 




Temperature [ K ] 



FIG. 5. Temperature dependence of the total entropy as 
calculated by thermodynamic integration from the ideal gas 
reference state. Only points in the temperature range where 
MD measurements have been performed arc shown. The ref- 
erence temperature T* = 380 K is also shown (dashed line). 




10'/T[K"'] 

FIG. 6. Energies of the inherent structures at the different 
densities as a function of 1 /T. The straight lines confirm the 
validity of Eq. (30) in the entire temperature range consid- 
ered. 
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FIG. 7. Anharmonic contributions to the energies, at the 
two indicated densities, together with the appropriate cubic 
fit Eq. (31). This contribution is integrated to directly calcu- 
late the anharmonic contribution to the vibrational entropy. 
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FIG. 8. (a) Density of states at fixed density pA at the tliree 
indicated temperatures. This quantity is the histogram of the 
square root of the eigenvalues of the Hessian calculated for the 
inherent structures, (b) Density dependence of the density of 
state at fixed temperature T = 320 K. The dashed line indi- 
cates the isosbestic frequency w* w 44 cm~^ at which all the 
curves intersect. The relevance of this feature is discussed in 
the text. 
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FIG. 9. (a) Temperature dependence of the average basin 
curvatures log{(jJk/(jJo)', this quantity, being a 

sum of logarithms, is very sensitive to the spectrum tails. 
Wo = 1 cm~^ sets the frequency scale, (b) Relation between 
the energy of the inherent structures and the average basin 
curvatures. The straight lines confirm the correlation between 
shape and depth of the inherent structures accessed by the 
system. 
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FIG. 10. Main panel: Harmonic contribution to the vi- 
brational entropy as calculated from the eigenvalues of the 
Hessian for the inherent structures. Inset: Anharmonic con- 
tribution to the vibrational entropy as calculated by integra- 
tion of the anharmonic contribution to the potential energy, 
as discussed in the text. 
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FIG. 11. Volume and temperature dependence of the con- 

figurational entropy Sc calculated as the difference between 
the total and the vibrational entropy. Solid lines are interpo- 
lations of the calculated points to Eq. (34). 




250 



300 



450 



500 



350 400 
Temperature [ K ] 

FIG. 12. Temperature dependence of the different contri- 
butions to the total entropy (closed triangles) at the fixed 
selected density pA'- harmonic (open squares), configurational 
entropy (closed circles) and anharmonic (open diamonds). 
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FIG. 13. Diffusion constants together with the correspond- 
ing power law fits (solid lines) predicted by the MOT. The 
breakdown of this prediction and the crossover to an acti- 
vated dynamics is evident. See text for a discussion of this 
point, (a) As a function of temperature, (b) As a function 
of the inverse temperature in order to stress the exponential 
dependence at the lowest temperatures. 
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FIG. 14. MOT paxameters as calculated from the diffusion 
constants. Main panel: Critical temperature Tc{V) (open 
circles) together with the value calculated in Ref. [42] (closed 
circle). The dashed line is only a guide for the eye. Inset: 
Power law exponent f{V). 




FIG. 15. Test of the 

Adam-Gibbs relation log _D(T) oc (l/TSc) for five different 
densities. Note that this linear relation holds both above and 
below the estimated critical temperatures Tc- 
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